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ABSTRACT 


The seasonal variation of permanent currents within 
Baffin Bay is investigated. A model based on physical reason- 
ing is developed, whereby the speed of surface currents may be 
calculated during winter months, when observations are not 
available due to ice coverage. Verification is made wher possi- 
ble utilizing late fall, and early spring data. [Limitations, 
applications, and variability of results are discussed. Two 
important applications are: first, to add to the oceanographic 
knowledge of the Baffin Bay area, and second, to provide an 
additional input to computer programs already in existence for 


forecasting the drift of sea ice within Baffin Bay. 
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1. Introduction 

Numerous oceanographic cruises have been made into Baffin Bay 
in recent years primarily by Canadian and United States oceano- 
graphic vessels. The purposes of these cruises are to extend 
man's knowledge of this region, and to explore its commercial 
and military potential. This research, and other traffic as well, 
is hindered in large part by the yearly cycle of winter ice growth 
which generally leaves Baffin Bay accessible to shipping only dur- 
ing the short summer season, July through September. Thus, while 
a great deal is known about Baffin Bay during these summer months, 
the winter oceanographic regime is almost totally unexplored. It 
is, then, the first intent of this paper, to consider one aspect of 
this winter regime, that is, the surface circulation pattern. 

One important aspect of extending knowledge of this area is 
the problem of forecasting the drift of sea ice. The purpose and 
usefulness of such forecasts is well established; they are present- 
ly being made in the Baffin Bay area on a routine basis by hand 
calculation. Both the Fleet Weather Facility, Monterey, California, 
and the U. S. Naval Oceanographic Office, are in the process of 
developing computer programs to accomplish this task; and, more 
specifically, Knodle (1) in 1964 developed a computer program for 
forecasting the wind drift of sea ice. Verification of test data 
utilizing this program indicated that,if the additional input of 
permanent currents could be introduced into the program, the accur- 


acy of results might be significantly improved. This input of 





permanent currents would be in the form of vectors at grid 
points of any scale desired and would be added vectorially to 
the present output, which considers only wind effects. While 
this final step of programming was not completed in this paper, 
the basic data for its completion are developed and examined in 
some detail, together with a discussion of the variability of 


the currents and the limitations of the methods employed. 





2. Background 

The first extensive attempt to determine the circulation pattern 
in Baffin Bay was made in 1928$ and,while much has been done since 
that time, the work of Kiilerich (1939) (2) on the Danish Godthaab 
expedition of 1928, still stands as the most authoritative account 
of the summertime surface currents that exist within the Bay. The 
results of this work, as modified by the more recent observations, 
is shown in figure 1. Exhaustive descriptive accounts of this 
region are available from any number of sources (3, 4)$ so only a 
few distinctive features which are pertinent to the problem at hand 
will be mentioned. Again, with reference to figure 1, it should be 
noted that Baffin Bay is largely an isolated body of water with 
limited exchange taking place through the narrow channels to the 
north, and across the Davis Strait to the south. A ridge there 
forms a definite sill, with a limiting depth of less than 730m. 
The continental shelf is narrow on both the Greenland and Canadian 
sides of the bay, with steep slopes falling off rapidly into the one 
large basin, whose depths are in excess of 2000m. 

The circulation is cyclonic with one small gyral to the immedi- 
ate north of Davis Strait and one large gyral centered over the 
Baffin Basin. Mean velocities are shown in cm/sec. 

Various climatological atlases show that,from mid-November 
through late May, Baffin Bay is virtually ice-bounds: and thus in the 


following section a model is presented to show how the surface current, 


at any given location within Baffin Bay, is diminished with the advent 


of total ice coverage. 
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3. Wintertime Regime 

Due to the fact that no winter measurements have been made, 
this section is devoted to the development of a model, which, when 
applied, will describe the wintertime variation of the surface 
currents within Baffin Bay. 

Two aSSumptions will be made initially: first, that the 
Summertime circulation is primarily maintained by the wind, and 
further, that the contribution made to the surface circulation from 
the Artic Basin and the North Atlantic is at a maximum during the 
summer, decreasing to a minimum by mid-winter. The latter part of 
this first assumption will be discussed in section 6. The former, 
is justified qualitatively on the basis of climatological wind data, 
and descriptive accounts of the area. For example, mean sea pres- 
sure charts for the summer months show Baffin Bay to be under the 
influence of a cold low pressure center. Further, the Polar Basin 
is under the influence of a cold high which contributes to steady 
easterly winds across the northern portion of Baffin Bay throughout 
the summer. The second assumption is that, once Baffin Bay becomes 
completely ice covered, no wind stress is transmitted through the 
ice. In other words, the contribution the wind makes towards main- 
taining the current is zero. Naturally, there will be some shifting 
of the ice due to tides, freezing effects, etc.; but in the mean the 


ice will be, for all practical purposes, stationary. 





If the assumptions above hold, the only real forces acting 
to alter the basic currents are first, the force of friction be- 
tween ice and water, and second, the frictional forces between 
the sea floor and the water. The latter is assumed to be rela- 
tively small and is not considered in this model. This is reason- 
able when one considers that the continental shelf is narrow, and 
that several investigators, including Kiilerich, found the deep- 
basin water to be near motionless. Further, if one assumes that 
the energy gained from currents entering the bay, is approximately 
equal to bottom and boundary losses, then the omission of these 
losses is further justified, at least qualitatively. 

The fundamental winter relationship can now be seen: the 
energy of the currents is diminished at the rate at which energy 
is dissipated by frictional stress between ice and water. If this 
rate is known, together with the vertical distribution of velocity, 
it follows that the rate at which the current velocity decreases 
can be calculated. In order to find this rate, first let the stress 
be given by 


Yak st (1) 


coefficient of friction between ice and water, 


where k. 


surface current velocity, defined to be the current 


u 
Found a few meters beneath the sea surface, and 


{ = density of the water at the same level 





The rate of work per unit area is Cy, which is to be 


equated to a rate of energy dissipation per unit area: 


x, ot 7am 
Cu = “ye (Kinetic Energy) oy. 


Area 
1 & lew: blige?” Prt Sar i bang 2urs 2) 
Tay = Ry, Me uo -~AG[e oly + le OOchALYTSU, QA ht (2) 


where as mean current speed in the friction layer, 


u, = the deviation from the mean U, 


Qi.= the mean density in the friction layer, 
and, ah = depth of the frictional layer. 
The corresponding parameters in the layer from the depth of frict- 
ion, to the level of no motion are Ua Q2 ,and Q&2 as shown in 


nena i 
figure 2. Note that uo = +<uU D is applicable in either 


layer. 
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Schematic representation of a 
wind-driven velocity profile 
showing quantities represented 
in equation 2. 


FIGURE 2 





From the observations made of velocity profiles shown in 
Appendix A, two initial observations can be made. First, the 
profile can change radically within a short period of time, and 
within short distances; and, second, these changes take place 
throughout the entire column of water. While these profiles may 
contain many inaccuracies as will be discussed in section 6, the 
observations do suggest that a loss or gain of energy at the sur- 
face influences the entire column of water that is in motion, 
Thus, the assumption is made that the kinetic energy loss is dis- 
tributed from the surface to the level of no motion nearly uni- 
formly. This means that the total loss, that is, the right hand 
side of equation (2), may be expressed solely in terms of the 


loss within the friction layer, or more specifically 


Rk. § & = 2 (ae + <u >Je ah) (3) 


where the constant @ is dependent upon the level of no 

motion (LNM) and Ah - The level of no motion applicable to 
Baffin Bay is discussed in section 4; but, if for example, the 
LNM is found at 1000m, and Sh = 50m, then C would equal twenty 
ete. (or C = LNM/ ah ). This doubtless is an oversimplification. 
This assumption, however, appears to fit the observations, and 


gives realistic results. 


Thus, in equation (3), the total kinetic energy loss is ex- 
pressed solely in terms of the loss that occurs within the fric- 


tion layer. This layer will be examined in some detail in the 





remainder of this section. If several basic relationships, taken 
from Shuleikin (5), relating to wind drift of sea ice are intro- 
duced, it will be shown that the thickness of the friction layer, 
re may be expressed purely as a function of the surface current 
velocity, and further that this same expression may be used for the 
model under consideration. Thus, for a current flowing beneath 


stationary ice: 


(a) sh= J Also 


) 


_ Tah 
@) us | AL WAR , 
(Cc) u = AVN on ' 
ee 
@) T= ka dV 
where: 
T = tangential stress between air and ice, 


aS 
i 


eddy viscosity coefficient, 


S 
l 


density of the air, 


a 
| 


= w sing, where w is the angular velocity 
of the earth's rotation, and = latitude, 


Hee coefficient of friction between air and 


water 
A = the wind factor, 
V = the wind velocity, and 


Ah , Sor UU, are as previously defined . 
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Tah sin? 


From (b) and (c) is obtained A= ia ie BV 


which when substituted into (a) yields 


L _ tT Asunp- _ tH Ra So.) {Send 
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Finally, solving (c) for V and substituting into the above, 


ah =Nuw, where (V= = bs. parc? - 


With sufficient accuracy kK. = 2x10 , and on the basis of studies 
by Several authors, A = 1.27x10°. Substituting these, and other 
known values, into the above expression yields 

| ah = 473u (4) 

where ah is in meters if u is in meters/sec. 

In the above development the tangential stress T operating 
on the sea surface was caused by wind. Shuleikin concluded, how- 
ever, that with ice coverage the frictional layer would also be 
present, but would be generated by the surface friction between ice 
and water. Further, he concludes that the numerical value of the 
coefficient N will not vary perceptibly because "it is known that 
the turbulent regime affecting the value sh becomes established 
under the influence of any velocity of the surface current u, 
independent of whether this velocity is generated by the tangential 
stress of wind origin, or by the stress caused by the friction of 
iceson the water", Thus, in this model equation (4) will be used to 


relate u and Ah 


Jal 





With the magnitude of the currents that we are dealing with» 
it is easily seen that this depth of friction will rarely exceed 
fifty or sixty meters$ and within this thin upper layer it can be 
assumed with a good degree of accuracy .that Sa =P , where Sg is 
the density of the water a few meters beneath the surface, and gi 
is the mean density within the friction layer. 

The numerical value of the coefficient of friction, Ko has 
been investigated in some detail by several authors, and in part- 
icular by Brown and Crary (6), Fukutomi (7), and Shuleikin (5). 
The values varied with location and with the age of the ice. 
Wittman and MacDowell () modified Shulekin's coefficients for use 
specifically in Baffin Bay, and Knodle (1) concluded after further 
investigation, that a mean value of .013 was the best approximation. 
This value of Ry then will be used for this model. 

Now, from equation (4) and the density approximation above ¢ 


equation (3) takes the form, 


: f mee ii an uA pi “ 4 \ 
Jan ao cal 2 eae a Hf) (5) 


It will be shown in the following paragraphs that the expression 
fa +<u>] may be treated, with some approximation, solely 
as a function of the surface velocity, u, if it is assumed that the 
Shape of the velocity profile, once established, does not alter 
significantly. While the wind-driven profiles given in appendix A 


show a great variety of shape, it should be noted that the greatest 





variability occurs when the currents are weak and in shallow 
water. Profiles associated with stronger currents in deep water 
(i.e., from S to 10 cm/sec. at the surface) show a definite 
uniformity of shape. It is argued here that most of the profile 
variability comes about as the result of non-uniform influences, 
such as wind gradients at the surface, summer run-off, inflowing 
currents, tides, or internal waves, and that these disturbing in- 
fluences are either totally absent. or at a minimum during the 
winter season. 

Considering the above, we will now construct a model for the 
current velocity profile. First, it is assumed that during the 
winter, the frictional component of the velocity (Ve) is analogous 


to that given by Ekman; then 


(ec) Ve(#) = Heh : 


where, ugh velocity found at the depth of friction, and 


\. (2) 


frictional component of the velocity found at 

depth 2 within the friction layer. 
Punthes it is assumed that,due to the proximity of land within 
Baffin Bay, there is no rotation of this frictional component with 


depth. It follows that: 


beth > Ve(2) 


where u, = the magnitude of the ¢urrent found at depth & 
within the friction layer. If, for example, at the time the ice 


becomes stationary, the depth of friction is found at 50m, then, by 


3 
equation (4), U ah may be calculated and a velocity profile plotted 
(see figure 3). 
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Several comments are needed concerning this profile. First, 
the velocity at 32 = ah is the proper veloaity to give the rela- 
tion between u and ah , i.e., ah = ae Ua . Second, equations 
(e) and (f) are being utilized as the winter profile model since no 
direct and detailed profile studies for flow beneath stationary ice 
were available to the author. Finallly, the introduction of this 
profile necessitates redefining u, the surface current velocity, 
as will become apparent in the following paragraphs. 

To find the relation between the surface velocity u, and ut in 
the friction layer, several velocity profiles were constructed utiliz- 
ing equations (e), (f), and (4) for assumed values of u,} , from 
12 cm/sec to 4 em/sec. For ease of computation each profile was 
approximated by a straight line as shown in figure 3. The velocity 
gradient with depth could then be treated as a constant (du/dz = C), 
and 

(z) “= Ue +c of 

(h) Us =4e + az 


where uc is the velocity at depth 5 meters, defined as the 
Surface current velocity. 


By definition 


(iy) uU = Ue % » 

sail u' = 2 (2-46) ) 

co oh 2 
Sage a ea 


and hence 


(3) 





By utilizing equations (g) and (j) to calculate the quantities 

Lu > and u for each velocity profile approximation, two important 
conclusions were reachedt first, that the quantity Cotecu' >] 
could be approximated quite accurately by fie ss <ul >]s Lt 7 
and secondly, that the same expression can be described in terms of 
the surface current velocity u, again without significant loss of 


accuracy 
si if ® , er 2 
[a+< u>]=1.10 = Kk, Us ; (Sa) 


where the constant kj; equals ¢.3¢. If equation (5) is modified 


in accordance with equation (2a), then 


kus = -£($ Us ]H73 Uh) , (6) 


With the rapid decrease in current velocity near the ice surface 
associated with the winter profile of figure 3, it becomes appar- 
ent that u. no longer adequately defines the velocity within the 
friction layer, but rather that Ue ah is the representative current 
velocity. Utilizing the same winter profiles as in the above 

deve lopment eae. it was found that the approximation Up, = 3.18u, 
was quite accurate, providing the range of velocity profiles 
chosen is not too large. That is, within the range utilized, 4-12 


em/sec for u,, , the accuracy is quite good, and would be the 


range of values expected for the major portion of Baffin Bay, Hence, 


if we define u = ugh , and substitute the above approximation into 


equation 6, then its solution 
Bu (4 -bo) 


u (4) = ulte) @ (7) 
where 
Ke = He ¢ (ae )ek 
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It is seen then, that at the time the ice becomes stationary the 
current velocity u of equation (/) represents the current found at 
the depth of friction, and that this depth will decrease as the 
winter season progresses approaching the surface by the time of 
spring ice break-up (see figure W) . Thus, u at the depth of frict- 
ion is applicable only during the time the ice is stationary, while 
u several meters beneath the sea surface (i.e., at 5m) is applicable 
for the remainder of the year. 

Thus, equation (7), shows the time variation of the current at 
depth ah , where u(tg) represents the current at the initial time 
of ice coverage and is known if the assumption is made that the 
current velocities of figure 1 are representative of ugh, at the 
time the ice becomes stationary. Hence, we may calculate the current 
at any time u(t) between initial ice coverage; and the time of brea&- 
up of the ice. The application of this equation is discussed in the 
following section, and a curve showing the seasonal variation of the 
surface current is developed utilizing this equation. 

The transition from the summer to the winter profile, is be- 
lieved to take place within a short time after the ice becomes sta- 
tionary. This is analogous, for example, to tidal currents near 
hottom, as given by Defant (9). In this transition period me Kew E . 
loss is confined almost entirely to the friction layer, and the re- 
duction of K. E. with depth does not occur until after the winter 


profile becomes established. This transition period is shown schemat 


ically in figure 4, together with a profile, say several months later. 


7 
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4. Application 


In this section the application of equation (7) developed in 
the previous section will be discussed, together with a qualita- 
tive treatment of the spring and summer current build-up, thus 
completing the mean yearly cycle. 

Two facets of equation (7/7) must be dealt with prior to its 
application. First, the time the ice becomes stationary and the 
time break-up of the ice occurs must be determined; and, second, 
the value of C,discussed briefly in section 3, needs to be calcul- 
ated. The latter will be treated first. Recall that C = LNM/ah 
and that the assumption has been made that kinetic energy is dissi- 
pated from the surface to the level of no motion nearly uniformly. 
This implies that iepievall of no motion descends with time. By 
equation (4) it is seen that the depth of friction ah is a 
function of the surface current velocity u and thus is likewise a 
function of time (see figure 4). Hence, a profile study was under- 
taken, one purpose of which was to investigate the variation of 
these quantities (see appendix A). From it, some approximations 
could be made. From the fall data, it was estimated that the 
level of no motion in deep water occurred at approximately 1000m 
and,as previously discussed, ah = 50m. Only a very few spring 
Soundings were available to the author, and only one in deep water. 
The velocity profile constructed from them is as shown in figure 5, 


together with a fall velocity profile taken at approximately the 
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same location. While this indeed is sketchy information, both 
the spring and the fall profiles,in this instance, lend support 
to earlier assumptions, and in particular to the assumption that 
an energy loss is dissipated with depth nearly uniformly. From 
this profile it was estimated that at the time of break-up of the 
ice LNM = 500m and Ah = 20m. Thus, C varies from the fall to 
the spring only from 20 to 25; a mean value of 22.5 is chosen for 
C in equation (7). 

The span of time during which this model would be applicable 
was determined on the basis of climatological data to be from 
October through April. The limitations of utilizing these dates 
are discussed in section 5; however, they should be approximately 
correct for an average year. Now all the quantities of equation 
(7) are known, and u(t) may be calculated utilizing the velocities 

the 
of figure 1 as u(t,). The expression from equation (7), SRL ( 4-49) 
represents the fraction of the summer time current velocity and 
may be plotted versus time with t, = 0 on 1 October and t (in 
seconds) up to 1 May. The plot is shown by the solid line portim 
of figure 6. 
The assumed comine and summer build-up of the current is 


shown by the dashed portion of figure 6. This build-up was 
handled in a purely qualitative fashion due to time limitations. 
Now, figures 1 and 6 together permit determination of the mean 


current velocity for any time during the year or location within 


Baffin Bay. 
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5. Limitations 
There are three major limitations to the winter-time model 
developed in the preceding sections which require discussion. 
First, the assumption was made, based on climatological data, 
that the ice cover became approximately stationary on 1 October, 
and that break-up occurred on 1 May. This statement implies that 
“these times occur Simultaneously throughout the Bay, when in fact 
the freezing and break-up occur over a period of time with vary- 
ing locations. An attempt was made to divide Baffin Bay in four 
approximately equal areas and develop separate curves similar to 
figure 6 for each quadrant. This trial was abandoned,however, as 
it introduced more uncertainties and complications into the pro- 
blem, and these disadvantages seemed to outweigh whatever benefits 
might be derived. For example, a declining current in one quadrant 
would certainly effect the next quadrant into which the current was 
flowing. No ready solution to this problem is proposed; it is felt 
that whatever modification is introduced would lead to only small 
changes in figure 6. The use of climatological data also implies 
that any one year is like the next, and this obviously is not the 
case. Figure 6 could, however, be quite easily reconstructed, 
based on information from local reconnaissance flights for any 
given year. 
The second major limitation concerns the variation in the 


Summer time build-up of currents. Evidence indicates that this 
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variation in some years might be large. For example,in 1964 
residual ice was observed in Melville Bugt throughout the summer, 
because the mean pressure patterns normally associated with this 
season were displaced enough from their normal positions that the 
winds, which normally would advect ice out of Baffin Bay and con- 
tribute to the current build-up, were simply not present. During 
Such an anomolous year, the model would obviously not work. How- 
ever, this unuSual situation is a rarity. 

The final limitation concerns the sum of the approximations 
and assumptions made in the development of the model itself. 
While these assumptions and approximations have all been discussed 
individually and justified, the cumulative effects could be addi- 
tive and result in serious error. Particularly in the assumption 
of the shape of winter velocity profile and in the treatment of 
C=LNM/ ah could large errors arise. However, in the absence of 
winter data, it is not possible to assign a number to this cumula- 
tive error. In section 7, verification of a forecast utilizing 
this model is quite good. This in itself is not proof that all the 
approximations and assumptions are valid, but it does at least lend 


Support. 
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6. Current Profile Study 

In order to verify the assumptions made in earlier sections, 
as well as to verify figures 1 and 6, an attempt was made to ob- 
tain the published results from the latest falls and earliest 
spring oceanographic cruises that had been made into Baffin Bay. 
These data were received from the Bedford Institute of Oceano- 
graphy, Darthmouth, N.S., via the Canadian Oceanographic Data 
Center (10, 11). Utilizing the standard techniques of computing 
and summing dynamic height anomolies, a number of velocity profiles 
was constructed. These profiles, together with their locations, 
are shown in appendix A. While, in large part, this study failed 
to yield the verification that was hoped for, some interesting ob- 
servations can be made. The remainder of this section is devoted 
to a discussion of what was hoped to be gained by this study, 
whether success or failure was met, and whenever possible why. 

First, it was hoped that the shape of the velocity profiles, 
at least for the fall data, would fit the preconceived notion of 
a wind-driven velocity profile (see figure 2). While in a number 
of instances, particularly in deep water, this was the case, a 
great variety of shapes resulted. While causes for variability 
were mentioned briefly in section 3 (tides, internal waves, etc.), 
it is believed by the author that a major limitation of the dynamic 
method of computing currents is that oceanographic station data can- 


not be taken simultaneously (that is, with only one ship). For the 


mS) 





profiles shown in appendix A, a lapse time of from two to six hours 
occurred between adjacent stations; if reorientation of the dynamic 
topography were to take place during this interval, distortion in 
the magnitude of the currant at a given depth would occur. 

It was hoped, secondly, that the level of no motion could be 
firmly established for the spring and the fall, as various investi- 
gators, and particularly Kiilerich, had stated that this level will 
vary from season to season, perhaps even irregularly from place to 
place. While some substantiation of the figures used in this model 
were gained from this study (see figure 5), it can be noted from 
appendix A. that this level cannot be fixed with confidence. Here 
also, one must be careful to distinguish between the reference 
level and the level of no motion. For this study, the reference 
level was chosen in every instance at the greatest depth that the 
adjacent soundings would permit. However, in some instances the 
soundings did not approach the bottom; and, hence, the profiles give 
only a velocity that is relative to whatever current exists at this 
reference level. This of course means also, that the true level of 
no motion could be at quite a different depth from that implied by 
the profiles. 

One of the initial assumptions made was that exchange of water 
between Baffin Bay and the bodies of water to the North and South was 
at a maximum during the summer, decreasing to a minimum by mid-winter. 


As can be seen from appendix A, only the exchange across Davis Strait 
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could be examined. While the assumption made cannot be validated 
on the basis of this study, it is interesting to note that, if the 
assumed level of no motion is correct, outflow greatly exceeded in- 
flow in both seasons with no inflow occurring with the fall data, 
and very little with the spring data. From previous calculations 
of current flow across Davis Strait, Dunbar (4) had stated that 
the volume transport iS approximately in the ratio of 2 to l with 
outflow being twice as great. He concluded therefore, that Baffin 
Bay received its waters in approximately equal quantities from the 
north and the south, but failed to consider runoff, which on the 
basis of this study would appear to make an important contribution. 
Hence, while the ratio of 2 to 1 would appear to be an underestima- 
tion, if, in fact, water does enter Baffin Bay in approximately 
equal quantities from the North and South, one could draw the tenta- 
tive conclusion that the amount of water entering Baffin Bay is 
small during all seasons. 

On the basis of the fall data, an effort was made to validate 
the current velocities given in figure 1. Of the twenty-five fall 
velocity profiles plotted in appendix A, fourteen, or 56%,were in 


good agreement with figure 1, that is,within + 3cm/sec at the sur- 


} 
face. One further limitation of the dynamic anomoly technique is 
that it gives only that component of the current which is normal to 


a ling drawn between adjacent stations. This surely accounts in 


large part for the variation from figure 1 in several of the profiles. 
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The remainder, however, and some are the exact opposite of what 
one would expect, can only be explained in terms of the limita- 
tions of the methods employed outlined in previous paragraphs. 
In conclusion, there was only fair agreement between the obser- 
vations and figure 1, with no consistent error being observed 
that would indicate that it should be changed. 

A verification of figure 6 was hoped to be accomplished on 
the basis of the spring data. These data, however, due to ice 
coverage were restricted primarily to Labrador waters, with only 
eight soundings taken within Baffin Bay. Of these eight, all 
but one showed a definite decrease in the velocity as the model 
predicts. The verification presented in the next section is an- 


other test of the model, and of figure 6. 
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7. Veritieation 

In Section 1, it was stated that one of the primary reasons 
for investigating the seasonal variation of the permanent cur- 
rents within Baffin Bay was that the results of the study could 
be used as an additional input to Knodle's computer program for 
forecasting the wind drift of ice, developed specifically for 
Baffin Bay. In this section, the results of Knodle's three veri- 
fication runs are shown in figures 7, 8, and 9, together with a 
modified forecast which includes the effects of permanent currents. 
This latter forecast was arrived at by adding to the wind-drift 
calculations of Knodle's, hand calculated, advection of the ice 
edge by currents in accordance with figures 1 and 6. While this 
is a laborious and time consuming task, the results were gratify- 
ing; and the job would be a simple one for the computer. 

As can be seen, there was indeed a significant improvement 
over Knodle's results in each of the three cases. One might con- 
clude, however, that the magnitude of the currents was too small, 
or in effect, that figure 6 overestimates the winter seasonal de- 
cline of the current. It is possible that melting might account 
for a large part of the remaining error. It is concluded on the 
basis of the foregoing that, while figures ] and 6 may be altered 
in the light of later findings, their use in forecasting the drift 


of sea ice within Baffin Bay tends to improve the drift forecasts. 
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8. Conclusion 

In the preceding sections, a model has been developed which 
predicts the winter decline of the permanent currents within 
Baffin Bay. In the absence of winter data, it was necessary to 
make a number of assumptions, on the basis of physical reasoning, 
about the nature of the flow. The spring and summer build-up was 
handled qualitatively thus completing the yearly cycle. A limited 
verification of the results was obtained. 

Tt is concluded that, despite the limitations of the model, 
steps have been taken toward adding to the oceanographic knowledge 
of the Baffin Bay area and providing the basic material for a com- 
puter input to forecasting the component of sea ice associated 
with permanent currents. Finally, it is recommended that the re- 
sults of this paper be combined with the results of Knodle. It is 
felt that there would result a fully operational program, far 


faster and more accurate than is presently done by hand calculation. 
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APPENDIX A 


VELOCITY PROFILES 


On the following pages are shown the velocity profiles as 
constructed from the data contained in (10) and (11), together 
with the locations for the fall profiles (figure A. Sept. l/7 - 
24, 1960) and the spring profiles (figure B. June 1 to ll, 1963). 

The computations were made in accordance with the Sandstrom- 
Holland-Hansen method of computing currents. It is assumed that 
the reader has some familiarity with the technique, so that only 
a nie Gocleee will be required here. 

l. Dynamic height anomolies were computed for the standard 
levels and summed upwards from the CS eee level (R.L.) 
to the surface in accordance with: AW 5 d Ap 
with the specific volume anomoly s , The R.L. was 
picked in all cases at the maximum depth that adjacent 
soundings would allow. 

2. Relative currents for each level were computed by the 


formula: 


7,4 Ae “aD ) 
UW), ~ 2 Scene (Ada ~ ie 
where 

W)= current found at a given level 


Wy. = current found at the reference level 
assumed equal to zero 


Ocean Oc0ecan 


Sto..A Sta..15 
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The results of these computations, as stated, are shown in 
the following pages together with other pertinent data such as 


the depth of the water, the time the sounding was taken, etc. 
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VELOCITY PROFILES AT STATION 
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VELOCITY PROFILES AT STATION 
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